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^ , Abstract 



We derive the form of the asymptotic series, as t — > +co, for a general 
solution h{t) of the non-linear differential equation h^t)"^ {h" {t)+h' (t)) = 
1. 



1 Introduction 

The purpose of this article is to describe the asymptotics of solutions of the 
second-order ordinary non-linear differential equation 

hit)%h"{t) + h\t)) = 1 (1) 

with initial conditions 

h{to) = ho>0, h'{to) = h. (2) 
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Before formulating the result let us describe our motivation and the origin 
of the problem which has its roots in the physical Hall effect. 

In a classical mechanics description the issue is to study the dynamics of 
a point mass moving in a periodic planar potential and driven by an exterior 
electromagnetic field where the magnetic field is constant and the electric 
field circular and created by a linearly time dependent flux tube through 
the origin, see for the origin of the model. The equations of motions are 
Hamiltonian. The time dependent Hamiltonian is 

H{t, q,p) = ^ip- A{q, t)f + V{q, t) on \ x 

with 

Here h and e are real parameters and V a smooth periodic function. In 
Newtonian form the equations of motion are 

q = E{q) + bBq- VV{q) in \ {0} 

where D is rotation by 7r/2 and E{q) = \q\-^]D)q= -dtA{q). 

We shall prove elsewhere that if b and e are nonzero the solutions are 
diffusive with or without direction depending on the direction of the fields. 
In this article we discuss the particular case when e = l,b = and = 0. 
In polar coordinates the Hamiltonian reads 

H{t, r, 0, pr, P4>) = \{p^r + -^ (P0 + 
and the equations of motion become 

P4> = 0, 0' = Pr = ^^^il^, r' =Pr. 

Consequently, p^ is a constant and r" = r"^{p^ + 1)^. After a shift in time 
we arrive at the equation 



The substitution 



r{t) =thi\nt) 
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leads to equation (|I]). 

In order to formulate our result we have to introduce some auxiliary 
notation. Let Sm,k G M[ai, 02, . . . ajt] be polynomials defined as follows: 

Sm,k{(ll, 0,2, • • • ; Cffc) = 0'nQ'i2 ■ ■ ■ (^imi 

ii+i2+---+im,=k 

m = 1,2, k = 0,1,2, ... . Clearly, Sm,k{o,i, 0.2, ■■ ■ , au) = if < m, and 
we set by definition so,a; = 5o,k- The polynomials obey the recursive rule 

fc-i 

Sm+i,fc(ai, 02, • • • , flfc) = ^ Smj(ai, 02, . . . , aj) at-j for m + 1 < /c. 

i=tn 

In the space of formal power series, it holds 

'00 \ ™ 00 

akX^ = ^ s„,fc(ai,a2, . . . ,afc)a;'', m = 0,1,2,.... 



X;afcX^/ = X;^a;^eR[[x]] 



^fc=l / k=m 

This implies that if 

00 

a = 

k=l k=0 

then 

00 00 k 

^ fma"" = '^9kX^ where fiffc = ^ /mSm,A:(ai,a2, ... ,afc) . (3) 

m=0 k=0 m=0 

Furthermore, set 

cr°(ai, 02, . . . , Ofc) = ^ ^ Sj^k{ai, ^2, • • • , ctfc) (4) 

j=i 

and 

cr^(ai,a2, . . . ,afc) = ^ f . J s^- ^(ai, 02, . . . , a/,) for m > 1. (5) 

Then it holds, in 

(00 \ 00 

1 + ^akx''] = ^al{ai,a2,...,ak)x'', 
k=l / k=l 

00 \ -™ 00 

l + ^a^a;''! = 1 + ^ o-™(ai, 02, . . . , a^) x'' for m > 1. 



fc=i / fc=i 
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Let {(3n}^=Q be a sequence of real numbers defined recursively, 

\ / ^ I — n o 1. C—n 



(6) 



j,k=0 j,k,e=o 

j+k=n+l j+k+e=n+l 



Here are several first values: 

R - 3 ^ _ 21 165 7245 

/i^i--^, /^2---, /^^3--^, /^4--^, .... 

For a fixed constant c e R we introduce a sequence of polynomials, 
Pn{c', z) e ^[z], n e Z+, by the recursive rule 

Po(c; z) ^2,z- c (7) 

and 

^k+l ^ 

Pn = 3(7°(po,Pl, . . . ,Pn-l) + 0-n-fc(PO,Pl, . . . ,Pn-fe-l) 

" (8) 

n 

This can be rewritten with the aid of the polynomials Sm,k-i 

Pn = : Sj^n{P0,Pl,- ■ ■ ,Pn-l) 

+ 2^ — ;^^rfc — )«j-,fe(po,Pi,... 

A;=0 i=0 \ J / 

For n > 1, the degree of Pn{c\z) is less or equal to n (this can be easily 
shown by induction when using the fact that for any monomial a^^di^ ■ ■ ■ a// 
occurring in a^{ai,a2, ■ ■ ■ ,ak) it holds '^ijSj — k). Here are several first 
polynomials Pn{z), 

Pi{c; z) = 9z — 21 — 3c, 

27 3 

P2{c;z) = -— ^2 + (90 + 9c)^- 228- 30c--c^ 

P3{c;z) = 27z^ - (^^ + 27c^ z2 + (1638 + 207c + 9c2)z 

69 

-3540 - 546 c - c\ 

2 

Now we are able to formulate the result. 
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Theorem 1. For any initial data (to, ^o? ^i) £ ]Rx]0,oo[xM there exists 
a unique solution h{t) to the problem jj^, ^ on the real line. Moreover, 
there exists a constant c = c{to, Hq, hi) G M such that it holds, for all n G Z+ 
and t +00: 

h(t) = (Atr (l + ± ^0{{^^y')\ (9) 

where 

m=l ^ 

The degree of qk{c; z) is less than or equal to k. 
Remarks, (i) Several first polynomials qk{c; z) are 

3 1 
qiic; z) = — z c, 

/ , 27 2 9 9 \ 21 3 3 2 

qoic; z) = z + 1 c\z c c, 

' ^ 512 V64 256 J 64 64 512 ' 

, , 189 3 

q-iic; z) = z — 

' ' 8192 



/ 135 189 \ 2 f 549 45 63 ^\ 

\ c] + { \ c H 

V1024 8192 J \102A 512 8192 J 



57 183 15 2 7 o 
■ c - — — - — — c. 



64 1024 1024 8192 
(ii) In the final step of the proof, in Subsection ^4.5| , we shall show the 
following invariance property of the asymptotic expansion. Set 



A^{c-t) = {4tY/Ul + J2 



ln(4t)) j 



-00. 



t^ 

k=l 

with n G Z+ and c, t G M. Then for all s G M it holds true that 
An{c-t + s) = An{c-As-t) + t^l^0(^(^^^ ^ ) ast- 

The remainder of the paper contains all necessary steps to prove Theo- 
rem 0. We shall proceed as follows. In Section ^ we show the completeness 
and derive the first term of the asymptotic series. In Section ^ we make 
use of the fact that the second order differential equation can be reduced to 
a first order differential equation and we investigate the asymptotic prop- 
erties of the latter equation. These results are used in Section § to derive 
the asymptotic properties of the original second order differential equation 
and this way we complete the proof of Theorem |l]. Section |^ contains an 
additional remark on the asymptotics of the Lambert function. 
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2 Basic properties of the differential equation 

The differential equation (|I]) is equivalent to the dynamical system 

{x',y')= (^y,^~y^ on M = ]0, oo[ xM. (10) 

Proposition 2. The flow of ^Tdj ) is complete and so for all initial data 
{to, ho, hi) G Mx ]0, oo[ xM there exists a unique globally defined positive so- 
lution h of ^ with initial conditions h{to) = ho, h'(to) = hi. 

Proof. We use the following criterion (c.f. |jl|, chap. 2.1.20]): The flow of 
a vector field ^ on a manifold M is complete if there is a proper map 
f G C^(M, R) which meets the estimate 

3A>0,B>0,WpeM, \^-f{p)\<A\f{p)\ + B. 

In our case M = ]0, oo[ xM, ^ = yd^ + {x^^ — y) dy and we choose /(x, y) = 
+ + x~^. With this choice we have 

|e-/(x,i/)| = \2xy-2y^\<x^ + 3y^<3fix,y). 

Moreover, for any bounded set 5* C M the inverse image f~^{S) is bounded 
and separated from the border of the half-plane M: there exists e > such 
that f~^{S) C [e,oo] X M. This implies that / is in fact a proper map and 
the proposition is proven. □ 

Proposition 3. Let {x{t),y{t)) , with t G [to!Oo[, be a solution of the dy- 
namical system 1^1^). Then there exists T G [to)Oo[ such that 

yt,s,t>s>T, V2{t + c(s))^/^ < x(t) < y{s) + V2{t + c(s))^/^ 

where c(s) = \ x{sy — s. 

Proof. Set (in this proof) g{x,y) = x~^ — y and 

giix,y) = -3^- ^ + y = - + g{x, y)-3 — 
X* x"* V X/ X 

on M. Thus x' = y,y' = g{x,y) and j-^g{x{t),y{t)) = gi{x{t),y{t)). 

We shall show first that y'{t) is negative for all sufficiently large t. Note 
that 

V(x, y) G M, g{x, y) < ^ gi{x, y) < 0. 
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Thus if y'{s) < then y'{t) < for all t > s. Hence it is sufficient to show- 
that there is at least one t such that y'{t) < 0. Suppose the contrary. Since 
it holds 

there exists s > to such that y{s) > 0. Then both x(t) and y{t) arc increasing 
positive functions on the interval [s, oo[ and, in addition, x{t) < y(t)~^f^. 
So the function g{t) = g{x{t),y{t)) obeys 

Vt>.. ,'W = -(l + 3||).«)-3f!<-3.M'/'<0 

which clearly contradicts the assumption g(t) = y'{t) > for all t. 

Let now T > to be such that y'{t) = x(t)^^ — y{t) < for all t > T and 
fix s > T. For any t > T we have {\x{t)'^y = x{t)^y{t) > 1. Consequently, if 
t > s then 

x{t) >V2{t+ -x{s)^ - sj . 

To show the other inequality set, for t > s, z{t) = x{t) — \/2(t + c)"^/^ 
where c — c{s). We find that 

{e'z'y = e\z' + z") = e* - (v^ {t + cf")-' + ^ + c)"^/^ j 



It follows that 



z'it) < e'-'z'is) + ^ e-* /* (li + c)"^/^ dii 
16 Js 



and 

3^2 '•^ 



z{t) < z{s) + (1 - e^-*) z'{s) + ^ ^ e" (i^ + c)'^'^ du dr 

< z{s) + z'{s) + ^ [\l-e^-') {u + c)-'/Uu 

16 Js 

< + + ^l(. + c)-3/^ 
= z{s)+x'{s). 

But z{s) = and so x{s) < y{s) + (t + c)^/^ □ 
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Corollary 4. Ifh(t) is a solution of (j^) on [to,oo[ with the initial conditions 
h{tQ) = hQ > 0, h'ito) = hi, then there exists T > t^ such that h'{t) > for 
allt > T. 

Corollary 5. Ifh{t) is a solution of (j^) on [tQ,oo[ with the initial conditions 
h{to) = ho > 0, h'ito) = hi, then 

= V2t^/^ + 0(1) ast^+oo. 

Remark. This means that if we restrict ourselves in what follows to the initial 
condition (^ with /ii > we don't loose the generality as far as the asymp- 
totics is concerned. Furthermore, owing to the invariance of the differential 
equation in time we can set to = 0. This fact will be used in the course of 
the proof. First we verify Theorem |l] for the particular case when to = and 
hi > and then, in Subsection we shall extend the result to the general 
initial condition. 



3 A reduced differential equation of first or- 
der 

In accordance with the remark at the end of Section ^ we assume that to = 
and hi > 0. The second-order differential equation equation (|1]) is invariant 
in t and this is why it can be reduced to a first-order differential equation. 
Actually, using the substitution h{t) = (G~"^(4t))^^^, zq = 4:/hQ and go = 
h^hi where 

G{x) 



ds 

hi 



'A 9 (!) 

we arrive at a first-order nonlinear differential equation, namely 

l-lzg{z)-z'g'iz)\giz) = l, g{zo) = go- (11) 



We shall carry out the computations relating (|I]), (^ and (|TT|) later, in Sub 



section [4.1| . Here we concentrate on the study of (|TT]) on the interval [0, zo 



assuming that zo> Q and go> 
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3.1 Domain of the left maximal solution g{z) 

We shall need two equivalent forms of the differential equation, namely 



1 / 1 \ 3o(z) , , 

g\z) = — 1 - — - - 12 
^ ^ ^ g{z)J A z ^ ^ 



and 

■■'^W)' = (l - ^) . (13) 
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Remark. In what follows we use repeatedly the following elementary argu- 
ment: if ip and (p are two differentiable functions on ]a,b[ and the equality 
tp{z) = (p{z) implies ip'^z) < ip'{z), for all z E]a,b[, then the two functions 
coincide in at most one point z G]a, 6[. 

Proposition 6. Let g{z) be the left maximal solution of pl|/ Then g{z) is 
defined and positive on ]0, zq], and 

lim^(2) = l. (14) 

Proof. Let 7 be the minimal non- negative number such that g{z) is defined 
on ]7,zo]- Equation ([TT|) clearly excludes the possibility that g{z) = for 
some z G ]7, 2:9 ]. So g{z) is positive on ]7, zq]. Our goal is to show that 7 = 
and (p^ holds true. We split the proof into six claims. 
(z) 3^e]^,Zo] s.t. g{0 > 1. 

Suppose that g{z) < 1, Wz e]7,2;o]. Then, by ([T^), g'{z) < on ]7, Zq] 
and so there exists limzi-y g{z) = gi with g^ < gi < 1. Hence, by minimality 



of 7, it should hold 7 = 0. According to (p!3D, 

.3/4 



{z"'g{z))' < =^ g{z) > go (^) , G ]0, zo 



:) 

a contradiction. In the following choose ^ G ]7, Zq ] to be the largest number 
such g{i) > 1. 

(ii) g{z) > 1, Vz G]7,e[. 

Actually, g{z) = 1 implies g'{z) < 0. 

(Hi) 'y = 0. 

If 7 > then, by (|13D, \{z^/^g{z))'\ < 7"^/^ < 00 for all z g]7,^[. This 
means that z^^'^g{z) is absolutely continuous on ]7, ^[ and so liuizi^ g{z) exists 
and is finite, a contradiction with the minimality of 7. 

(zv)3r]e]0,^[s.t. g'iz)^0,\/ze]0,v[. 
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Equalling the RHS of (|T2[) to zero we get a quadratic equation with respect 
to g{z). Its solution is a couple of functions, 

2 



1 + 1- v/l-32' 

defined on the interval ]0,^[n]0, |[. Clearly, ip'i{z) > and ^2(2) < ev- 
erywhere on that interval. This implies that the RHS of ([T^ ) vanishes in a 
point 2; from that interval if and only if g{z) = ifi{z) or g{z) = (f2{z) and 
in such a case either = g'{z) < f'i{z) or = g'{z) > ^2(2). Thus fii^z) 
coincides with g{z) in at most one point z, and the same is true for (p2{z). 
Consequently, there exists r] G]0,^[n]0, |[ such that g'{z) doesn't vanish on 
]0,?7[. Choose r] having this property. 
(v) g'iz)>0,\/ze]0,v[- 

If g'{z) < 0, \fz ^]0,ri[, then gi = lim^^o exists (finite or infinite) and 
91 > 9{v) > 1- Oil the other hand, in virtue of (|12|), 



z'g\z) = l--^^-^zgiz)<0. (15) 

Equality (|13]) implies that (z^^'^g^z))' > on ]0,r][ and so 

/^N 3/4 

giz) < giri) i^j) , Vz e]0,v[. 
Consequently, limzio z g{z) = 0. Sending 2; to in (^) gives 

1 - — < 0, 

a contradiction. 

(vi) lim^|o5'(^) = 1- 

From Claim (v) follows that gi = limzio g{z) exists and 1 < (71 < g{i])- 
Suppose that gi > 1. Then one concludes from (|T2|) that there exists 6 G ]0, ri[ 
s.t. 

g'{z) > 4, Vz G ]0, 6[ , where d=]-(l- — ] > 0. 
^ ^ \ 9i J 

This implies 

a contradiction. □ 
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Corollary 7. The maximal solution g{z) satisfies the integral identity 

g{zf = 2z''/^ r s-'/\g{s) - 1) ds, \fz G ]0, ]• (16) 
Jo 

Proof. Rewrite (|12D as 

(.^/vy = 2.-^/^(^7-1), 

and integrate from to z when taking into account (0). □ 

3.2 Asymptotics of the left maximal solution g{z) at 
z = 

Lemma 8. If P > then 

(z^) as z 10. 



J z 



Proof. 



1 1 ' 



< e ^ ( — e ) + i2z)^e ^ ( — 



□ 



Notation. (a)o = 1, («)j = a{a + 1) . . . (a + j — 1). 
Lemma 9. Let 5 > 0, z/ G M. Then the expression 

z-^e-^^' / s^-^e^^'ds. 



regarded as a function in the variable z, has the asymptotic series, as z I 0, 

oo 

Proof. It suffices to show that it holds, for all n G Z^., 

^5 rt— 1 n— 1 



/Q It J_ It J. 

j=0 j=0 

+ (z/)„e-^ / s'^+^-^e^ c/s . 

J z 
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Actually, according to Lemma M this means that the relation 



/5 n-l 



holds true for all n > —u, and so for all n E N. 

We shall proceed by induction in n. For n = this is a trivial equality. 
The induction step n n + 1: 



5 rS 



j-5 



□ 

Let us consider equation ([Tl|) (without the initial condition) in the space of 
formal power series C[[2;]]. Its solution 

oo 

~g{z)=Y,a,z' eCM] (17) 

fc=0 

is unique, with the coefficients being determined by the recursive relation 

ao = 1, ak+i = ( 2^ + 4 ) for /c > . (18) 

Several first coefficients are 

3 15 483 

ao = 1, «! = -, 02 = y, as = ■■■ ■ 

Proposition 10. The left maximal solution g{z) has an asymptotic series, 
as z I 0, that is equal to 

oo 
k=0 

Proof. We have to show that, for all n G 

n 

g{z) = Y,(^kz' + o{z^). (19) 

k=0 
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We shall proceed by induction in n. The case = means that \im.ziog{z) 



1 and is covered by Proposition |g. Let us suppose that (|T9D is valid for some 
n G Z+. Denote (in this proof) 

a_=liminf^^^l^%l^. 

Similarly, a+ designates the limes superior, as 2 | 0, of the same function. 
Thus the induction step n — > n + 1 means to verify that 
We shall do it in three steps. 

(i) < a-a+i < a+. 

If 6 < a_ then there exists 5 > Os.t. g{z) > ELo ^-^ e]0,5[. 

Further, from the recurrent relation ([l8|) one finds that 

00 / k \ 00 2 

fc=0 V j=0 / fc=0 2 

and the assumption (^) implies that the asymptotics of g{zY is given by a 
truncation of the formal power series giz)"^, namely 

n 2 
fc=0 2 

Combining this with ([l6| ) leads to the conclusion that, V2; G ]0, 

Hence b < a„+i for all b < a- and consequently a_ < ttn+i- The inequality 
a+ > ttn+i can be proven symmetrically. 

(a) Let A > be a fixed parameter. In the second step we introduce an 
auxiliary function f{z) whose choice depends on whether n = or n > 0. In 
the former case we set 

i 3 1 / 5 1 



(p{z) = —z 4e Az y s 4eAs ds 



31/ 5 1 



(Xz) 4e Az / s 'ie^ ds (21) 
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and in the latter one 

.1 n+l 



^V^fcS^+'^-^gl (22) 

- k=0 



where 



V^fc = Ofc - ( ^ - 7 ) «fc-i for 1 < < ra, 



4 ; - ^ ^ (23) 

Observe that = ai — |a;o = 0. In the case n = 0, {p{z) solves the 
differential equation 

^'{z) = ^Mz)-l)-l^, (24) 
Xz^ A z 

while in the case n > 0, ^{z) solves 

^' W = -M^) - 1) - - E V'.-^-^ • (25) 

z^ 4 z ^-^ 

fc=2 

We claim that, in the both cases, the asymptotic behaviour of V9(z), as ^ | 0, 
is given by 

n 

ip{z) = J2 o^kz^ + A + o(^"+') . (26) 

A:=0 

Actually, equality ( [26| ) follows directly from Lemma y. In more detail. Lemma 
^ gives, for the case n > 0, 



1 — n I— n V / 1 



(/?(Z) = 

fe=o i=o 

n+l m 



jn=0 fc=0 ^ ^ '^-'^ 

The coefficients yj^, as given in (p3l), have been chosen so that the asymptotics 
( ]26| ) is satisfied. This is to say that equalling (^) to (|26| ) leads to a system of 
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linear equations on the coefficients ipk whose unique solution is exactly (p3| ) 
as follows from the identity 

J,— Q \ / m—k k=l ' ^ ^ m—k 

The case = is even more straightforward. 
(iiij a_|_ < ttn+i and a_ > a^+i- 

Let us show the first inequality, the other one can be proven analogously. 
We shall need the asymptotics of the function 



' \k=0 J \k=l 

n-1 

k=0 

Again, assumption (0) implies that Ylk=o ^^^^^ ^ truncation of the power 
series 

-/ X 2 = {9-^)9+-; 

g[z) z^ 4 z 

= -(q^)'-q' + -- -. 

Here we have used that g solves (p!2D . Combining (|l^ and (|20|) one arrives 
at the formula 

lk = ak+2- [k + ^Aak+i. (28) 



Now we can compare the functions g{z) and <^{z). Let us choose A > 1 in 
(|21| ) resp. (|22|). Suppose that ip{z) = g{z) at some point z. Then, owing to 
and (H, it holds 

^'(z) - g\z) = Y^Mz) - l){g{z) - A), 
Az g[z) 



when n = 0, and using (^Sf ), (|23D and 



. n+l 

^'i^)-9'iz) = -^Mz)-lf-Y.^kz'-' 

' k=2 

n—l I „^ ^n—l^ 



(l-A)a„+iz"-i + o(^"-i) 
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lim inf — = a_ - A a„+i < (1 - X)a„+i < 0. 



when n > 0. In any case, there exists 6 > s.t. ^p{z) = g{z) imphes 
(p'{z) < g'{z), Wz g]0,5[ (in the case n = 0, we need also that g{z) > 1 if ^ 
is sufficiently close to 0, see Proposition |^). This means that the functions 
g{z) and (p{z) coincide in at most one point z G ]0, 6[. Furthermore, we have 
already shown that a_ < and so, using also (|26|) , we conclude that 

gjz) - f{z) 

Thus there exists a sequence s.t. z„ i and g{zn) < ^p{zn), Vn. Con- 
sequently, g{z) < ip{z) on a right neighbourhood of 0. Hence, in virtue of 

a+ s nm sup r— — A 

The claim is a consequence of the limit A | 1. □ 

Corollary 11. The left maximal solution g{z), after having been defined at 
z = by g{0) = 1, belongs to C°°([0,zo])- 



Proof. Observe that consecutive differentiation of equation (|T2D jointly with 
Proposition 10 imply that, for any m E z^~^^ g^''^\z) has an asymptotic 
series at z = which we shall call YlT=o have to show that g G C", 

Vm, and this in turn amounts to showing that = for k < m + 1. Let us 
proceed by induction in m. The case m = was the content of Proposition 
^. Assume now that g G C"". Then = for k < m + 1, and the mean 
value theorem implies that 

lim inf ^('"+^) (z) < '^3^'"]^^+^ = < sup ^(™+^) (z) . (29) 

zlo dz zlO 

On the other hand, since z"^^"^ g^'^~^^\z) has an asymptotic series the limit 
liuizio g^"^~^^\z) always exists and equals either ±oo or 0:^+2 depending on 
whether there exists an index k < m + 2 s.t. o;^^^ 7^ or not. However the 
property (p9|) clearly excludes the ffist possibility. □ 

4 Asymptotics of a solution h{t) of the second 
order differential equation 

Except of the last subsection, we still consider the particular case when to = 
and hi > (see the remark at the end of Section We shall proceed to 
the case of general initial condition only at the very end of the proof, in 
Subsection ETS, 
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4.1 Reduction of the second order differential equation 



Let us now complete some computations concerning the reduction of the 
second order differential equation (|l]), (||) to a first order differential equation. 
Let g{z) be the left maximal solution of the first order differential equation 

1 - ^zg{z) - z'^g'{z)^ g{z) = 1, g{zo) = go, 

where 

4 

^0 



2^0 — T4) 9o — h^hi 

fin 



^From Section |^ we know that g{z) is a positive function from the class 
C°°([ 0,2:0]) (Proposition 1^ and Corollary |ri|). Consider the function 



G{x) = I -jj^ , h'^<x<oo. (30) 



ds 
9 (!) 

Then G G C°°([ /iq, oo[ ), G is strictly increasing, Giji^) = 0, and, owing to 
([T^), \imx^ooG{x) = 00. So the inverse function satisfies G'^ G C°°{[0, oo[) 
with G-\Q) = hi. Set 



hit) = (G-i (At)) 



1/4 



Then h(t) solves the problem (|I]), (||). 

Actually, G(/i(t)^) = At, G'{h^) = g{A/h^)-^, and so 



'^''^'^^^^«(i|. (31) 



4G"(/i4) dt " \h^ 
Differentiating (|3l|) once more gives 

3h\h'Y + h^'h" = -16 h-'^h'g' (^-^ 

Denote for brevity z = A/h'^. Hence 

/i3(/i' + h") = g{z) - ?>h-\h'^h'f - 16 h-\h'^h')g\z) 

3 

= 9{z) - -zg{zf - z^g{z)g\z) 
= 1. 
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Furthermore, h{0) = {G-\0)f^^ = ho and 



h'{0) = G-\0)-'/' j-G-\0) = (^^G{ht 
4.2 Asymptotics of G{x) 

First let us find, in C[[2;]], the reciprocal element to the formal power series 
K^) = EZo'^kz" defined in (0), (^). Set 



n=0 



The formal power series g{z) solves the differential equation (|12D and so an 



easy calculation shows that f{z) = g{z) ^ solves the differential equation 

m = -\f{zni-f{z)) + li^. 

z^ 4 z 

On the other hand, this differential equation implies a recursive rule on the 
coefficients /3„, namely rule (^) preceding the formulation of Theorem |T]. 

Lemma 12. The asymptotic series at infinity of the function G{x) defined 



in (|5^ is given by 



where 



~x-3 ln(x) + c-4V^ f - 

k \x 

k=l ^ 



r ( ^ -l + ^]ds-ho' + 3 ln{h^). (32) 



Proof. It holds 



oo 



1 1 3 , ^ 
+ - \ ds. 
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According to Proposition |T0| we have the asymptotics at infinity, 

k 



1 , 3 /4V 



9 (!) ' k=2 

The claim then follows straightforwardly. □ 

4.3 Asymptotics of G^^{x) 

Let us now focus on the inverse function G^^ . 



Lemma 13. There exists such that for all x > x^, it holds true that 

G 

Proof. Choose > G^Hq) > so that, for all y > y 



0<G^\x)-x<—^(x-G(x)). (33) 
x — 4 



< 1 ^ <- and y- G{y) > 0. (34) 



This is possible owing to Proposition |T0| and Lemma |T2[ Set = max{4, y^}. 
For X > x^ fixed define a sequence {yn}^=o by the recursive rule 

yo = X, yn+i = x + yn - G{yn). 
Owing to our choice, yn^ x > x^, > y^, for all n. Using (^) one finds that 

rvn+i ryn+i / 1 \ , 

In virtue of (0), the sequence satisfies 

4 

< yn+2 - yn+1 < 4(ln(?/„+i) - ln(?/„)) < - {y^+i - yn)- 

X 

Since yi — yo = x — G{x) we get 

/4\" 

< yn+1 - yn< {-] {x- G{x)), \/n. 

By the choice of we have 4 < x and, consequently, the sequence {?/„} is 
convergent. The limit y = lim?/„ solves = x — G{y) and so y = G^^(x). 
Moreover, 

00 

< y-yn = ^{yk+i - yk) 

k=n 

k=n ' 

The particular case n = in this relation is nothing but our claim. □ 
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Combining Lemma |T3| with Lemma |T2| one immediately gets 

Corollary 14. It holds true that, as x ^ +00, 

G~\x) = x + 0(ln(a;)). 

Recall that in @, (||) we have introduced polynomials cr^(ai, 02, . . . , a^) 
labeled by indices m > and k>l. 

Proposition 15. For all n G Z+ it holds true that, as x 00, 

G-(x) = . + Mc; 1„(.)) + ± Pi<^ + 0( (!!^)"" ) (35) 

k=l ^ \ ^ J 

where the polynomials Pn{c', z) have been defined in and the constant 

c is given by equality ^^). 



Proof. Corollary |TJ implies 

'ln(x) 



,n(G-.w)^,„W.o(!l^M)._L_.iH-o(; 

Combining (^6|) with Lemma ^ one derives the relation 



(36) 



X 



(37) 

valid for all n > 0. Setting n = in ( pTl) one arrives at the case n = in 
(|35|) . To finish the proof one can proceed, in the obvious way, by induction 
in n when repeatedly using relation (|37|). □ 



4.4 Asymptotics of h{t) for particular initial data 

We already know that h{t) = {G~\4:t)f^^ solves the problem h{tf{h"{t) + 
h'(t)) = 1, h{0) = Hq, h'{0) = hi. Using the known asymptotics of G~^{x) 
we get 



and consequently 

hit) = iity" (1 + ± + o( (1^)" 



1/4 



k=l 
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where 



in=l ^ ^ 

So gfc are exactly the polynomials introduced in Theorem It is also easy 
to see that the degree of qk{c; z) is less than or equal to k since the same is 
true for the polynomials p„ with n > 1 and degpo = 1- This observation in 
fact completes the proof of Theorem |^ in the case when to = and hi > 0. 



4.5 General initial conditions 

Consider first a solution h{t) of with the initial conditions h{0) = ho, 
h'{0) = hi, assuming that hi is positive. Then, as we already know, the 
asymptotic behaviour of h{t) is described by Theorem |I], i.e. equality (^) 
holds true with c = c(0, ho, hi). Choose s G M and set h{t) = h{t + s). Then 
h{t) solves equation ([l|) and satisfies the initial conditions h{0) = ho = h{s), 
h'{0) = hi = h'{s). But h'{s) is positive for s sufficiently small and so 
equality (HD applies to h(t) as well, with c being replaced by c = c(0, ho, hi). 
Equating the asymptotics of h{t + s) to the asymptotics of h(t) one arrives 
at the equahty 



k=l 




(St 



valid for t +00 and every ri G Z+. /^From p8| ) it is not difficult to derive 
the relation between c and c, it reads 

c = c - 4s. (39) 

Thus the invariance of the differential equation (||) is reflected in an invariance 
of the asymptotic expansion of its solutions, as expressed by relations 



(|39|). It is also clear that these relations must hold true not only for s small 
but even for all s G M. 

Choose now arbitrary initial data {to, ho, hi) G Mx ]0, oo[ xM and let h{t) 
be the corresponding solution. Then, as we know from Corollary^, h'{t) > 
for all sufficiently large t. Fix s > to such that h'{s) > and set h{t) = 
h{t + s). We use once more the already proven fact that h{t) satisfies equality 
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d^), with c being replaced by c = c{0, h{s), h'{s)). This imphes that the 
asymptotic behaviour of h{t) = h{t — s) is given by 



n.+l 



n G Z+. But in that case one deduces from (|38|), (|39D that satisfies 
equahty @ as well, with c = c + 4s. Theorem |I] is proven. 

5 Additional remark: comparison with the 
asymptotics of —W-i{—e~^) 

This is a digression whose aim is to emphasize a rather close similarity of the 
asymptotic behaviour of the function G~^{x) with that of Lambert function. 
The Lambert function W{z) gives the principal solution for w in z = we^ 
and Wk{z) gives the k^^ solution. Surprisingly, it is not documented in some 
standard text books and reference books on special functions though we may 
have missed some sources. On the other hand, the Lambert function seems 
to have attracted even in a rather recent period some attention, particularly 
from the computational and combinatorial point of view (see for a sum- 
mary). It is also implemented in some standard computer algebra systems 
like Maple and Mathematica where it is called LambertW and ProductLog, 
respectively. Let us just briefly recall that W{z) is analytic in a neighbour- 
hood of z = with the convergence radius equal to e~^, 

k=l 

The coefficients have a combinatorial interpretation when counting distinct 
oriented trees. 

Consider now the equation 

y - My) = X, 

or, equivalently. 

It is elementary to see that for x g]1, +oo[ there are exactly two real solu- 
tions, yi{x) and 2/2(2;), with yi{x) €]0, 1[ and y2{x) g]1,+oo[. The both 
solutions can be expressed with the aid of the Lambert function, namely 

yiix) = -iy(-e-^), y2ix) = -Vr_,(-e-^) . 
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The aim of this remark is to point out that the asymptotics of the second 
solution, i.e. —W-i{—e~^), as x ^ +00, can be derived in a way quite similar 
to what we have done in Subsection [4.3| when treating the function G~^{x). 
To this end let us recursively define polynomials Pk{z), 

Po{z) = z, Pk+i{z) = a^^^{po{z),pi{z), ...,pk{z)). 

For A; > 1, the degree of the polynomial Pk{z) is k. Here are several first 
polynomials: 

Pl{z) = Z, P2{z) = Z-^Z^, p^{z) = Z - ^Z^ + ^Z^, .... 

Proposition 16. It holds, as x +00, 

-iy_i(-e-^) = a; + 0(lna;) 

and, for n > 0, 



-w'-.(-0- + E^ + o((!!i£) 



k=0 

The proposition can be proven using a similar approach as the one used 
in the proof of Lemma |13|. In fact, this asymptotic expansion is well known 
and is in agreement with what has been published in Q and though 
the derivation and presentation here is somewhat different. 



Acknowledgements. R.D.B. wishes to thank FONDECYT (Chile) 199- 
0427, and the action C94E10 of ECOS-CONICYT. P.S. wishes to gratefully 
acknowledge the partial support from Grant No. 201/01/01308 of Grant 
Agency of the Czech Republic. 



References 

[1] Abraham R., Marsden J.E.: Foundations of Mechanics. Addison- Wesley, 
1978. 

[2] Halperin B.: Phys. Rev. B 25 (1982) 2185ff. 

[3] Corless R.M., Gonnet G.H., Hare D.E.G., Jeffrey D.J., Knuth D.E.: On 
the Lambert W Function. Advances in Computational Mathematics 5 
(1996) 329-359. 

[4] de Bruijn N.G.: Asymptotic Methods in Analysis. North-Holland, 1961. 
[5] Comtet L.: C. R. Acad. Sc. Paris 270 (1970) 1085-1088. 



23 



